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We give general expressions for the vector asymmetry in the angular distribution of protons in 
the nonmesonic weak decay of polarized hypernuclei. From these we derive an explicit expression 
for the calculation of the asymmetry parameter, oa, which is applicable to the specific cases of iHe 
and ^aC described within the extreme shell model. In contrast to the approximate formula widely 
used in the literature, it includes the effects of three-body kinematics in the hnal states of the decay 
and correctly treats the contribution of transitions originating from single-proton states beyond the 
s-shell. This expression is then used for the corresponding numerical computation of oa within 
several one-meson-exchange models. Besides the strictly local approximation usually adopted for 
the transition potential, we also consider the addition of the first-order nonlocality terms. We find 
values for oa ranging from —0.62 to —0.24, in qualitative agreement with other theoretical estimates 
but in contradiction with some recent experimental determinations. 

PACS numbers: 21.80.-|-a, 13.75.Ev, 21.60.-n 
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I. INTRODUCTION 

While the A hyperon, in free space, decays 99.7% of the time through the mesonic mode, A —> ttN, inside nuclei 
this is Pauli-blocked, and, already for A > 5, the weak decay is gradually dominated by the nonmesonic channel, 
NN —> NN, where the large momentum transfers involved (« 400 MeV/c) put the two emitted nucleons above the 
Fermi surface. This decay mode is interesting since it offers a unique opportunity to probe the strangeness-changing 
weak interaction between hadrons. For a recent review of hypernuclear decay, see Ref. 0. 

For a long time, the experimental data for this process was restricted to the full nonmesonic decay rate, Tnm-, and, in 
some cases, also the partial ones, r„ = r(An ^ nn) and Fp = r(Ap —> np). More recently, the first results, obtained 
at KEK iili , for another important observable of nonmesonic decay, namely, the intrinsic asymmetry parameter, 
a A, are becoming available. This is experimentally more demanding, as it requires measuring the asymmetry in the 
angular distribution of protons emitted in the decay of polarized hypernuclei. On the theoretical side, however, oa 
carries important new information, since it is determined by the interference terms between the parity-conserving 
(PC) and the parity-violating (PV) proton-induced transitions to final states with different isospins. In opposition to 
that, the decay rates depend only on the square moduli of the separate components of the transition potential, being 
dominated by the PC ones. One expects, therefore, that the asymmetry parameter, besides being more sensitive 
to the PV amplitudes, will have more discriminating power to constrain the proposed mechanisms for nonmesonic 
hypernuclear decay. 

Most of the theoretical work on this decay mode constructs the transition potential by means of one-meson-exchange 
(OME) modeR the most complete ones including up to the whole ground pseudoscalar and vector meson octets (tt, r], 
K, p, Lo, K*) [a,|al3l- Recently we have extended such models to take into account the kinematical corrections due to 
the difference between the lambda and nucleon masses and the first-order nonlocality terms . There are also OME 
models that consider additional effects, such as correlated-two-pion exchange Q and direct-quark interaction |l(l| . 
In all these cases, to which we will refer below as strict OME models, the weak coupling constants for the pion are 
empirically determined from the free mesonic decay, and those of the remaining mesons by means of unitary-symmetry 
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arguments [1,0 • 

All such models reproduce quite easily the total nonmesonic decay rat^ ^nm = r„ + Fp, but seem to strongly 
underestimate the experimental values for the n/p branching ratio, r^/Fp |l|. However, there are recent indications, 
based on the intranuclear cascade model, that this might be due to contamination of the data by secondary nucleons 
unleashed by final state interactions (FSI) while the primary ones are traversing the residual nucleus El Another 
serious discrepancy between theory and experiment in nonmesonic decay concerns the asymmetry parameter. The 
measurements favor a negat ive value for and a positive value for ^He. However, all existing calculations based 
on strict OME models find values for oa between —0.73 and —0.19 0,^. Also, when results for 

are available in the same model, t hey are very similar to those for ^He, as expected, since the intrinsic asymmetry 
parameter, oa, has been defined in such a way as to subdue its dependence on the particular hypernucleus 
considered. A recent attempt ^ to explain this discrepancy along similar lines to those used for the n/p problem 
has failed. As might be expected, the FSI do have an effect in attenuating the asymmetry, but show no tendency to 
reverse its sign. The only theoretical calculations that attain some agreement with the experimental data for oa are 
a first application of effective field theory to nonmesonic decay El and a very recent extension of the direct-quark 
interaction model to include sigma-meson exchange |l li|| . However, in both cases, one or more coupling constants are 
specifically adjusted to reproduce the experimental value of oa for aHc. 

Most calculations of the asymmetry parameter make use of an approximate formula (Eq. (I43II . below) which, 
however, is valid only for s-shell hypernuclei. Since an essential aspect in the asymmetry puzzle presented above 
concerns the comparison of its values for ^He and it would be of great interest to have a simple expression 

that is applicable to both cases. This is the main objective of the present paper, in which a general formalism for 
the asymmetry parameter in nonmesonic decay is derived and subsequently particularized to these two hypernuclei. 
We start by presenting, in Section m the main steps in the derivation of the general expression, Eg. of the 
vector hypernuclear asymmetry in terms of decay strengths. This is equivalent to Eq. (27) of Ref. |13ll . However, 
we deviate considerably from that reference from this point onwards. The main difference is that we do not make 
use of spectroscopic factors, but rather rely on spectroscopic amplitudes, which can then be computed in the nuclear 
structure model of choice. This has, in our view, two great advantages. Firstly, the spectroscopic amplitudes can be 
determined without any ambiguity as to their phases. This is particularly important for the asymmetry parameter, 
where, differently from the case of the decay rates, one is dealing with an interference phenomenon. Secondly, since, 
due to the large value of the momentum transferred in the fundamental process, nonmesonic decay is not significantly 
affected by details of nuclear structure, one can choose to work in the extreme shell model. Doing this, much of the 
summation over the final states of the residual nucleus can be explicitly performed, leading to very simple expressions 
for the asymmetry. (See, for instance, Eq. 14211 .1 The scheme for computing the decay strengths by means of an 
integration over the available phase space is presented in Section llTTI and the summations needed for the asymmetry 
parameter are performed in Section ESI Finally, the numerical results obtained by applying this formalism to the 
calculation of ua for, both ^He, and in several strict OME models, are presented and discussed in Section m 
where we also summarize our main conclusions. Details of the derivation of the final expression for a a are given in 
Appendices EHO and some identities that have been used for this purpose are listed in Appendix IdI 


II. VECTOR HYPERNUCLEAR ASYMMETRY 


Single-A hypernuclei produced in a reaction, under favorable kinematical conditions, are known to end 

up with considerable vector polarization along the direction normal to the reaction plane, n = {Pt^+ x Pk+) l\PTr+ ^ 
Pk+\i of which they retain a significant amount, Py, even after they have cascaded down to their ground states by 
electromagnetic and strong processes [liEi Therefore, the initial mixed state from which the hypernucleus will 
decay weakly can be described by the density matrix US Eq.(9.29)] 


p{Ji) 


1 

‘IJj + 1 



3 

<// + 1 


PyJl ■ h 


( 1 ) 


where J/ is the hypernuclear spin. 

The angular distribution of protons emitted in the proton-induced nonmesonic decay, Ap —> np, of the pure initial 
hypernuclear state | J/M/) is given by Fermi’s golden rule as 


dF(J/M/ ^ p2tp) 


dPt 


P2 


= J dQ.p^ j dF ^ lipiSiU. 

S\S2Mp 


P2S2tp VpJpMp\V\JlMi)/ 


( 2 ) 


Here, PiSi and P 2 S 2 are the momenta and spin projections of the emitted neutron and proton, respectively, and we 
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have introduced the compact notation (Ti = c = \) 


/ 


dF... 



, |Pi +P 2 I 

2M 2M 2Mf 


-A 


Up Jp 


( 3 ) 


M being the nucleon mass; Mf, that of of the residual nucleus, which is left in state |ff Jf-^f) where vp specifies the 
remaining quantum numbers besides those related to the nuclear spin; and A^pj^, the liberated energy. (To avoid 
confusion, we will be using Roman font (M,m) for masses and italic font {M,m) for azimuthal quantum numbers.) 
Also indicated in Eq. 0 are the isospin projections = —1/2 and tp = +1/2 of the neutron and proton, respectively. 
The transition amplitude includes both the direct and the exchange contributions, f.e.. 


{piSitn P2S2tpVFJFMF\V\JlMi) 

= {piSitnP2S2tpiyFJFMF\V\JlMi) - {p2S2tpPlSitniyFJFMF\V\JlMi) , (4) 


where the round bras stand for simple (nonantisymmetrized) product states for the emitted nucleons and the transition 
potential, V, is extracted from the Feynman amplitude for the direct process 0. 

It is then possible to show ^ , by taking the appropriate average of Eq. that the angular distribution of 
protons from the decay of the polarized mixed state described by Eq. m has the form 


dT[p{Ji) P2tp 


dTt 


'P2 


(1 + PvAvP 2 ■ n) 

47r 


( 5 ) 


where Tp is the full proton-induced decay rate, and the vector hypernuclear asymmetry, Ay, is given by 


3 Em, MjajJjMi) 
+ 1 Em, 


( 6 ) 


The new quantities introduced above are the decay strengths, 


a{JiMi) 


J ditp^ J dF ^ \{PlSltnP2S2tpFFJFMF\V\JlMi)p,h.f. 


SIS2MP 


( 7 ) 


where the subscript p.h.f. indicates that one is dealing here with the transition amplitude in the proton helicity frame, 
in which the direction for angular momentum quantization is that of the proton momentum. Equivalently, one can 
choose, for the calculation of the decay strengths, a coordinate system having the z-axis pointing along the proton 
momentum, and proceed as usually. This is depicted in Fig.^ 

It is clear that, with the help of Eq. 0, one can extract the value of the product PyAy from the counting rates 
parallel and opposite to the polarization direction, by taking the ratio of their difference to their sum. Assuming 
that Py can also be independently measured, or calculated, this experimentally determines the vector hypernuclear 
asymmetry. Ay. 


III. DECAY STRENGTHS 


To compute the decay strengths, it is convenient to rewrite the transition amplitudes in Eq. 0 in the total spin 
{S, Ms) and isospin (T, Mp) basis. We start from the relation 

\P1S1tnP2S2tp) - \P2S2tpP1S1tn) 

= Y. (l/2si l/2s2|5Ms)(I/2t„ l/2tp|TMT) 

SMsTMt 

X [ \pPSM sTMt) - \-pPSM sTMt) ] , (8) 

where we have also changed the representation to relative and total momenta, given respectively by 


P = ^{P2-Pl), 
P = P1+P2- 


(9) 








4 



FIG. 1: Coordinate system for the calculation of decay strengths. 


Since we are taking = —1/2 and tp = +1/2, we can write 


(1/2 l/2tp\TMT) 



Sto) SmtO , 


and performing the transformation m in Eq. 0 , we get 


{PlSltn P2S2tpVFJFMF\V\JlMi) 

= - ^(1/2 Si 1/2 s2\SMs) {pPSMsTvfJfMf\V\JiMj) , 

SMs T 


( 10 ) 


( 11 ) 


where we have defined 


{pPSMgT ffJfMf\V\JiMi) 

= -^[{pPSMsTvfJfMf\V\JiMi) - {-pPSMsTvfJfMf\V\JiMi)] , (12) 

V2 

dropping, for simplicity, the Mt = 0 labels, as shall be done henceforth. Finally, introducing Eq. into Eq. 0 , 
and making use of the orthogonality of the Clebsch-Gordan coefficients in spin space, we are left with 


a{J[M[) = j dQp^ j dF 


SMsMf 


jpP^MsT vfJfMf\V\JiMi)p±,{, 


(13) 


For the integration in Eq. there are 6 momentum variables involved, namely, the components of pi and P 2 - 
These, however, are not all independent. The choice of z-axis in Fig. ^ two angular variables. Also, the 







5 


energy conservation condition in Eq. © gives one relation to be satisfied. This leaves 6 — 3 = 3 independent variables. 
A convenient choice is 


independent variables: p 2 , Op^ , 4>pi ■ 


Simple trigonometry, applied in Fig.^ leads to the relations 

= pI+pI- 2 piP 2 cos , 

P 1 +P 2 + 2piP2 cos , 




COS0p 

cos Op 


P 2 — Pi cos 0; 

2^ 

P 2 + Pi cos 0; 

p 


(14) 


(15) 


which, together with the energy conservation condition, determine all momentum variables in terms of the set in 
Eq. drill . Notice that the azimuthal angles of the several momenta are related as follows 


4>p — 4>pi + I (j>p — 4>pi ■ 


(16) 


For the transition amplitude, we expand the final state in terms of the relative and center-of-mass partial waves of 
the emitted nucleons 0, (2.5)], getting, 

{pPSMsT vfJfMf\V\JiMi)p±,{, 

= (47r)^ +'^)‘^YLiOp,(l)p,)]xp. 

ILXJ 

X {\plSMs\JMj){JMjJfMf\JiMi){pIPL\SJTvfJf]Ji\V\Ji), (17) 


where the values of p and Mj are fixed by the relations Mj = Mj + Mp = p + Mg + Mp- Due to the rotational 
invariance of V, the last matrix element in Eq. is independent of M/, and this label has, therefore, been omitted. 
For the same reason, the subscript p.h.f. has also been dropped. Notice that, from Eq. m and the well known 
behavior of the spherical harmonics under parity, one has 

{pIPLXSJTffJf; Ji\V\Ji) = ^ [l - (-)'+s+^] (j,IPLXSJTffJf; Ji\V\Ji) . (18) 

v2 

Upon integration on the angle (pp ^, Eq. m gives, then, 


ai JiMi) = J dcosOp, J dP ^ ^ (-)^ ^) ® 0)]a;. 


SMsMf ILXJT 


X {XpSMs\JMj){JMjJfMf\JiMi){pIPLXSJTffJf-, Ji\V\Ji) 


(19) 


It can be shown quite generally that 

[YiiOp, (j)p) 0 Yf{0p, (j)p)]xp [41'(0p, (jip) ® Yl>{9p, <(>p)]v^ 

= (47r)-i ii'LL'X' 

X ^ k{m'Q\kPi){LQL'Q\KQ){X' pkQ\Xp) 


UKk 

(i v k ] 

X {l L' K }[Yk{9p,ct>p)®YK{ep,^p)U, 

[X X' K j 


( 20 ) 


where I = \/2l + 1 and similarly for other angular momentum labels, 
performing the summations on Mg and Mp, Eq. m becomes 


a{JjMi) = ^(47r)'^ / dcosOp^ f ^ 

t/ 1 / 


T+T' 


Therefore, upon opening the square and 
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Y i^_Y+S+J+J'+Ji+Jf li'iPxX'JJ' 


ILXJ I'L'X'J' 


X Y KJi{JiMiKO\JiMi){lOl'{)\kQ){LQL'Q\K{))[Yk{ep,TT)(^YK{ep,Q)]^o 


UKk 


Jl K Jl 
J Jf J' 


K J' J 
S X X' 


I V k 
L L' K 
X X' K 

f r>T f \ f a jfrTir 


X (pIPLXSJTffJf-, Ji\V\Ji){pI'PL'X'SJ'T'ufJf-, Ji\V\Ji)* . 


( 21 ) 


IV. ASYMMETRY PARAMETER 

In order to carry out the summations on Mj needed in Eq. we first rewrite it in the form 

Ay = 3 

where we have introduced the decay moments 


Jl f^i{Ji) 
Jl + 1 (Jq{Ji) ’ 


CTo(J/) = 


Mi 


CTi(J/) = i~T Mia{JiMi). 

V Ji\Ji + 1) M, 


Mi 

Then we take advantage of the particular values 

(J/M/OOIJ/M/) = 1, 

(J/M/10| J/M/) = Ml/\J Ji{Jj + 1), 

and use the orthogonality relation 


Mi 


to get, for K = 0 and 1, 


\T+T' 


^k(.Ji) = -(47r)'‘j|K ^ / dcos9p, dF Y (“) 

J J grprp, 

X ^ ^ i-l'-P-l-L (_)A+S+J+J' + Jj + Jf ii'lPxX'JJ' 

ILXJ I'L'X'J' 

X ^ (Z0r0|fc0)(L0L'0|A0) [Yk{9p,F)®YK{ep,Q)U 


kK 


Jl K Jl ^ ( K J' J 
J Jf J' \ \ S X X' 


I I' k 
L V K 

K 

1' DT f \ f O jf'T^L 


X {pIPLXSJTvfJf]Ji\V\Ji){pI'PL'X'SJ'T'vfJf;Ji\V\Ji)* 
From 13, (2.13)] (see also My- 

{pIPLXSJTffJf] Ji\V\Ji) 

= {alal)^\\FFJF)MiplPLXSJT;jMp), 


ip 


( 22 ) 

(23) 

(24) 


(25) 

(26) 


where Ja = u-a l\ ja and jp = Up Ip jp are the single-particle states for the lambda and proton, respectively, and 

M{plPLXSJT-,iAip) = ^[l-{-y+SYiplPLXSJT\V\iAipJ). 


(27) 


(28) 


( 29 ) 
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We are working in the weak-coupling model (WCM), where the hyperon is assumed to stay in the Ja = lsi /2 single¬ 
particle state, and the initial hypernuclear state \Jj) is built by the simple coupling of this orbital to the core, taken 
as the ground state \ Jc), \Jj) = \{}aJc)Ji)- Under these circumstances, the two-particle spectroscopic 


amplitudes in Eq. lOSI) are cast as 

(J/ll = (-) 


J-\- Jl Jp 


JJi 


j" fp jf 


(30) 


To continue, we will adopt the extreme shell model (ESM) and restrict our attention to cases where the single¬ 
proton states are completely filled in \ Jc)- This is so, within the ESM, for the cores of, both ^He (J/ = 1/2, Jq = 0), 
and {Jl = l,Jc = 3/2). In such cases, the final nuclear states take the form li/pjp) = \{ip^Jc)JF), and we can 
associate the extra label vp with the occupied single-proton states, jp. Consequently, on one hand, only one term 
contributes to the sum in Eq. and, on the other, the corresponding single-proton spectroscopic amplitude in 

Eq. ® is given by 


{Jc\\4^\\FpJp) = {-Y- + ^^+^^Jp. 


(31) 


Notice also that, within this description, the liberated energies are independent of Jp, i.e., 
suggests rewriting Eq. as 


Aj^. This 


dF... = 


{2nf 


Jp 


Jp 


where 


J K is+ ii + 


f Pi , P2 , |Pl +P2P ,, 

' ^ ^ 2Mf 


(32) 


(33) 


Putting all this together and performing the summation on Jp, we finally get 


-XJi) = E 

ip STT' 

X ^ ^ i-i'-F-i-L ^_-^x+s+Ji+Jc-]p+>^ ii'LL'Xi'Yj'^ 

ILXJ I'L'X'J' 

X ^ (Z 0 r 0 |fc 0 )(L 0 L' 0 |K 0 ) [Yk{ep,Ti)®YK{9p,Q)U 


kK 


j\ Jl Jc 
Jl JA K 


K JA JA 

jp J J' 


K J' J 
S X X' 


I I' k 
L L' K 
X X' K 


X M{plPLXSJT-,]A}p)M*{pl'PL'X'SJ'T'-,iAip)- 

From Eqs. (EH) and m, it can be seen that all the dependence on J/ is contained in the factor 


3J/ 


1/2 Jl Jc 
Jl 1/2 1 


Ji + 1 j 1/2 Jl Jc 
Jl 1/2 0 


■// + 1 


for Ji = Jc + 1/2, 

for Ji = Jc-1/2. 


(34) 


(35) 


Thus, within the framework of the WCM, one frequently introduces the intrinsic A asymmetry parameter, oa, UM 
defined as 


OA = 


Av 

Jl + 1 A 
Jl 


for Jl = Jq + 1/2 , 

for Jl = Jq — 1/2 , 


(36) 


which does not depend on the hypernuclear spin, as we have just shown within the ESM and for core states having 
no open proton subshells. For such cases, we get 


OA = 


//I 

Wo 


(37) 









with, for K = 0 and 1, 


u>K = 8^2 ^ f dcos 0p^ ( dF-i^ X! 

jp •' STT' 

X ^ ^ i-i'-L'-i-L ^_^A+s+jp+i li'LL/xx'j^y 

ILXJ I'L'X'J' 


X 

X 

X 


Y, m'0\k0){L0L'0\K0) [Yk{9p,TT) ® YKi9p,0)U 

kK 


Uq 

J' 



v 

L' 

k 1 
K ^ 

}\s 

A 

A'/l 

U 

X' 



K 1/2 1/2 

jp J J' 

M{plPLXSJT-]!^p)M*{pl'PL'X'SJ'T'-\^\p). 


(38) 


It can be shown 0 that loq = Fp. Therefore, the new information carried by a\ comes from the numerator in 
Eq. 03, *-e-, from uji- 

The orbital angular momenta in Eq. 13HI obey the restrictions: 


^_y+l'-k ^ 

(-)L+L'-K ^ 

{_)k+K-^ = + 1 . ( 39 ) 


The first two follow from well known properties of Clebsch-Gordan coefficients with vanishing azimuthal quantum 
numbers. The last one can be obtained by first deriving the general relation 

[Yk{ep,(t>p)<S)YK{ep,ci)p)]i^ = (-)'=+'^-"+'^[n(0p,</p)®rK(0p,</>p)]._p. (40) 

Then, recalling that any spherical harmonic with azimuthal angle equal to, either 0, or tt, is real, one gets, 

[Yk{0p,Tr)®YK{9p,O)U = [Yk{0p,7r) (^YKi9p,0)Uo . (41) 


from which the third one of Eqs. (P|l follows immediately. The presence of the phase factor i ^ ^ ^ in Eq. 

may seem disquieting at first sight. However, by taking the complex conjugate of that equation, interchanging the 
dummy variables ILXJT ^ VL'X'J'T', and making use of Eqs. and of the symmetry properties of angular- 
momentum coupling and recoupling coefficients, one easily gets the relation uj* = uJk., showing that these quantities 
are real, as they should be by definition. 

To compute the two-body matrix elements defined in Eq. (123, we resort to a Moshinsky transformation [23 of 
the initial state, and phenomenologically add initial and final short-range correlations. (For more detail on this and 
related points, see Refs. Q and jj.) For ^He, the sole contribution to Eq. comes from the lsi /2 proton state, 
and one can put L = L' = K = 0. On the other hand, for also the lp 3/2 state contributes, in which case L and 
L' can each take the values 0 and 1. Consequently one could, in principle, have RT = 0,1 and 2 in Eq. But we 
prove in Appendix ^ that the contribution with K = 2 vanishes identically, both for k = 0, and for k = 1. Similarly, 
we prove in Appendix ini that the contribution with K = 1 vanishes for k = 0. We do not have an analytical proof 
that the contribution with K = n = 1 vanishes, but we show in Appendix [C] that it is, in any case, negligibly small. 
Therefore, only the term with K = 0 survives in Eq. (123, and it reduces to the following expression, that can be used 
for the two hypernuclei: 


= 



Y J d cos 9p^ J dFj^Yf,o{9p,0) 

Jp 


X 

X 


E E 

TT' LS IXJ I'X'J' 


^_ 'jX-\-X^ 


WxX'Tj'^ {101'0\k0) 


( K 1/2 l/2 \ ( K J' J \ ( l' I K 1 
{jp J F j{s X X' j{x X' L j 


X MiplPLXSJT;}Aip)M*ipl'PLX'SJ'T';}A}p), 


(42) 
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with L = 0 for the lsi /2 state, and L = 0 and 1 for the lp 3/2 state. 

It is interesting to observe that the presence of the Clebsch-Gordan coefficient in Eq. for K = 1, ensures that I 
and V have opposite parities. Since the initial state in the two matrix elements has a definite parity, this implies that 
all contributions to uji come from interference terms between the parity-conserving and the parity-violating parts of 
the transition potential. Furthermore, the antisymmetrization factor in Eq. ll^ shows that the two final states have 
T ^ T'. These are general properties of the asymmetry parameter, as mentioned in the introduction. 


V. NUMERICAL RESULTS AND CONCLUSIONS 


TABLE I: Results for the asymmetry parameter, aa, based on the nonmesonic decay of aHs- See text for detailed explanation. 


Model / Calculations 

<^o(lSl/2) 

t^l(lSl/2) 

fflA 

TT 

Strictly local 

Plus corrections 

0.5176 

0.6492 

-0.2254 

-0.2913 

-0.4354 (-0.4351) 
-0.4487 (-0.4456) 

{■k,V,K) 
Strictly local 

Plus corrections 

0.3322 

0.3920 

-0.1878 

-0.2412 

-0.5652 (-0.5852) 
-0.6153 (-0.6384) 

TT + p 

Strictly local 

Plus corrections 

0.5011 

0.5937 

-0.1227 

-0.1776 

-0.2449 (-0.2665) 
-0.2991 (-0.3155) 

{n,r],K) + {p,u},K*) 
Strictly local 

Plus corrections 

0.5352 

0.5526 

-0.2739 

-0.2974 

-0.5117 (-0.5131) 
-0.5382 (-0.5388) 

Experiment 

UA = Av(iHe) 

KEK-PS E278 

KEK-PS E462 [4] (preliminary) 


0.24 ±0.22 

0.11 ±0.08 ±0.04 


Shown in Tables Hand QH are the results obtained in the calculation of the asymmetry parameter, oa, based on the 
expressions of the previous section applied to ^He and ^a^, respectively. The values of ujq and uji are in units of the 
free A decay constant, tJ^^ = 2.50 x 10“® eV, and, in the case of ^aC, we give in separate columns the contributions of 
the lsi /2 and IP 3/2 proton states. Also included are the measured values for oa obtained from some recent experiments 
on the nonmesonic decay of these two hypernuclei. The value for oa in extracted from experiment KEK-PS E508 
and given in TableirTIwas taken from the preprint version of Ref. 0, since only a weighted average for and 
is explicitly reported in the published version, its value being —0.20 ± 0.26 ± 0.04. 

We consider several OME models, and for each one we give the results of two different calculations. First, are 
those corresponding to the strictly local approximation for the transition potential, usually adopted in the literature. 
Secondly, are the ones obtained when we add the corrections due to the kinematical effects related to the lambda- 
nucleon mass difference and the first-order nonlocality terms that we have discussed in Ref. Q. The first thing to 
notice is that these corrections act systematically in the direction of increasing the absolute values of all the tabulated 
quantities. The effect is typically in the range of 20-30% for Wk but only 5-10% for oa, tending to be larger in the 
t: p model and smaller in the complete model. To put it shortly, if one wants precise values for the asymmetry 

parameter, the correction terms should be included in the transition potential, but, in view of the present level of 
indeterminacy in the measurements, they can be dispensed with for the moment. 

In the case of aHo, we have also included, between parentheses, in Table Q the values for oa obtained with the 
approximate formula usually adopted in the literature Q, namely, 

23i[V^ae* - b {c* - V2d*) + Vi f {V2c* + d*)] 

|a|2 + |6|2 + 3(|c|2 + |d|2 + |e|2 + |/|2) 


Av(iHe) 


(43) 
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TABLE II: Results for the asymmetry parameter, ua, based on the nonmesonic decay of ^aC. See text for detailed explanation. 


Model / Calculations 

l^o(lSl/2) 

‘^o(1J33/2) 

0^l(lSl/2) 

t^l(lP3/2) 

OA 

TT 

Strictly local 

Plus corrections 

0.4111 

0.5206 

0.4724 

0.5954 

-0.1830 

-0.2400 

-0.1990 

-0.2596 

-0.4324 

-0.4477 

{■K,ri,K) 
Strictly local 

Plus corrections 

0.2788 

0.3336 

0.3161 

0.3817 

-0.1580 

-0.2057 

-0.1707 

-0.2217 

-0.5526 

-0.5975 

n + p 

Strictly local 

Plus corrections 

0.4138 

0.4922 

0.4607 

0.5514 

-0.0984 

-0.1461 

-0.1096 

-0.1596 

-0.2379 

-0.2930 

{■K,rj,K) + {p,uj,K*) 
Strictly local 

Plus corrections 

0.4391 

0.4619 

0.4803 

0.5083 

-0.2300 

-0.2546 

-0.2378 

-0.2599 

-0.5088 

-0.5303 

Experiment 

OA = -2Av(^iC) 

KEK-PS E160 [2] 

KEK-PS E508 \£ (preliminary) 



-0.9 ±0.3“ 
-0.44 ±0.32'’ 


“This result corresponds to an improved weighted average among several p-shell hypernuclei [J p.95]. 
^See text. 


where 


a = (np,^So|E|Ap,'^So), 

b = i(np,'^Po|E|Ap,'^So), 

c = (np,''*Si|P|Ap,'^Si), 


d = -{np,^I)i\V\Ap,^Si ), 

e = j(np,'^Pi|P|Ap,3Si), 

/ = -i{np,^Fi\V\Ap,^Si). 

(44) 


The extra factors in the transition amplitudes in Eqs. are due to differences in phase conventions, as explained 
in Ref. Q. It is important to emphasize that Eq. 63 is only an approximation, that can be adapted from the 
corresponding expression for the two-body reaction pn pA in free space |22j . As such, it ignores the fact that the 
final state of nonmesonic decay is a three-body one and the ensuing kinematical complications should be properly 
dealt with, which requires a direct integration over the available phase space as done in the expressions used here. 
More importantly, Eq. m does not include the full contribution of the transitions coming from proton states beyond 
the s-shell, being therefore of only limited validity, and should not be used for p-shell hypernuclei such as or, 
even worse, for heavier ones. This being said, comparison of the corresponding values for oa in TableOlshows that the 
formula works well within its range of validity. This conclusion is in agreement with our preliminary result reported 
elsewhere |23). which was restricted to one-pion-exchange only. 

Coming now to ^a^, it is evident in Table HII that the p-shell contributions to ujq and are by no means negligible, 
being in fact of the same order as those of the s-shell. However, they are also in approximately the same ratio, 
so that the effect on a a, given by Eq. m, is much smaller. This corroborates the theoretical expectation that 
the intrinsic asymmetry parameter, oa, should have only a moderate dependence on the particular hypernucleus 
considered. Presently we are investigating to which degree this remains true for more general cases, such as that of 
“B m. Notice that in the previous section we have explicitly proven that oa is independent of the hypernuclear 
spin. However this does not, by itself, exclude the possibility that it might depend on other aspects of hypernuclear 
structure. 

In closing, let us remark that we have derived simple formulas for the evaluation of the asymmetry parameter, 
which exactly include the effects of three-body kinematics in the final states of nonmesonic hypernuclear decay and 
correctly treat the contribution of transitions originating from proton states beyond the s-shell. As to our numerical 
results, let us first of all observe that the calculated values of oa in the four OME models considered here vary from 
—0.62 to —0.24. This broad spectrum of values indicates that the asymmetry parameter can indeed be a powerful tool 
to discriminate between different theoretical mechanisms for nonmesonic decay, requiring for this purpose, however, a 
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more precise experimental determination of this observable than those presently available. Secondly, the fact that, for 
each of these OME models, the results for ^He and are very similar is compatible with the general expectation 
that oa should depend little on the hypernucleus. Finally, the negative value systematically obtained for oa for the 
two hypernuclei indicates, once again, that it will be hard to get a positive or zero value for it in the first case, at 
least within strict OME models. The puzzle posed by the experimental results for a a in s- and p-shell hypernuclei 
remains unexplained. 


APPENDIX A-. K = 2 CONTRIBUTIONS TO wo AND wi 


As mentioned in Section im to compute the transition matrix elements M. appearing in Eq. (EHl), we perform a 
Moshinsky transformation of the initial Ap state 0 Eq.(A.l)], 

^ f 1/2 JA 1 

UaJp'.!) = jAjp ^ ^ ‘S' < Ip 1/2 jp / 

AS [ A 5 J J 

X (^nl NLX\nAlAnplpX) \nlNLXSJ) 

nlNL 

= C{nlNLXS;iAipJ) \nlNLXSJ) , (Al) 

nlNLXS 

where [nl NLX\nAlAnplpX) are the Moshinsky brackets with their phases adapted so as to conform with our 
convention for the relative coordinate as discussed in Appendix A of Ref. Q- We have put bars over I, X and S 
to distinguish them from the analogous angular momenta in the partial waves of the final NN state. This is not 
necessary for L and J, since all transitions are diagonal in these two quantum numbers. Introducing Eq. m into 
Eq. 123), one gets 


M {plPLXSJT; jAjp) = Y JAjp M {plPLXSJT; nl N XS ), 


nlN \S 


where 


M{plPLXSJT-,nlNXS) = ^ [l - (-)'+®+'^] (plPLA5'jr| U |nZ AL AS'J) . 

v2 


The transition potential can be decomposed as 




(A2) 

(A3) 

(A4) 


where the isospin factor A is equal to 1 or Ti • T 2 , for isoscalar or isovector interactions, respectively, and the fli are 
rotationally invariant operators having definite spin and spatial ranks, z.e., operators of the form 


n, = [A/*(cri,cr2)(g)S/‘(r,V)]oo . (A5) 

Due to the algebraic properties of the Pauli matrices, the rank can be at most equal to 2. For the OME models we 
consider, including eventual kinematical and nonlocality corrections, the several possibilities are: 


PC terms < 

\ 1^1 = 0 

for central and r ■ V forces 
(spin-independent or spin-spin), 


= 1 

for linear spin-orbit forces, 


1 

[ i^i = 2 

for tensor forces; 


PV terms 


for all kinds. 

(A6) 

As can be seen in Eqs. (A.3) to (A. 15) of Ref. Q, 

the different terms have matrix elements of the general form 


iplPLXSJT\viW^ \nlNLXSJ) 

= G,{ILXSJT;IXS) {PL\NL) {pl\vi{r)di{r)\nl ), 

(A7) 
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where {PL\NL) are the overlaps of the c.m. radial wave functions and the di{r) are, either unity, or one of the 
effective differential operators defined in Eq. (A. 16) of that reference. The important point is that the Gi are purely 
geometrical factors, involving, at most, 3j and symbols. Scrutinizing these equations more closely, one notices that 
the dependence on A, A and J can be isolated as follows 

G^ilLXSJT- rxs) = {-y I ^ ^ ^ J ^ 


This result is completely general and depends only on the application of the Wigner-Eckart theorem to operators of 
the form given in Eq. (IA5II . (See, for instance. Chapter 7 of Ref. [^, or Section lA-5 of Ref. 

Taking these ideas into account in Eq. (P|l. it is clear that the summation over A, J, X' and J' can be performed 
first, leading to a remaining summand that is proportional to 


EE( 

A,7 A',7' 


)A+J+J' 


X 


X 


X 


K 1/2 1/2 

Jp J J' 


K J' J 
S X X' 



' 1 

1' 

\ 

k 


L 

L' 

K 

J 

[a 

X' 

K 

/ 


[ Ia 1/2 JA 

\ Ip 1/2 jp 

[x S J 


Ai^iAMAi^iAl 
I L I j \s J S j 


( Ia 1/2 JA ' 

\ Ip 1/2 jp > 

[ A' S' J' ^ 


X' Vi' X' 
V L' I' 


X' vp X' 1 
S' J' S J ' 


(A9) 


Actually, the A is always in the lsi /2 state, and one can make use of Eq. (6.4.14) of Ref. to replace the above 
expression by 


X = 


X 


X 


X 


^(_)A+J+J' j2<^Sj2jS 


A,7 A'J' 


K 1/2 1/2 

Jp J J' 



S' 1/2 1/2 


Jp 


L J' 



(AlO) 


where we have dropped an irrelevant factor. 

For a proton in the s shell, the Moshinsky transformation requires that L = L' = 0, and the 9j symbol in Eq. (rnm 
selects AT = 0 as the only possibility. For a proton in the p shell, there are two alternatives for the relative and 
c.m. angular momenta, namely. 


p shell 


alternative 1: 1 = 1 and L = 0 , 
alternative 2: 1 = 0 and L = 1, 


(All) 


and similarly for the primed quantities. In principle, therefore, there are three possibilities for K, namely, AT = 0,1 
and 2. 

It is clear from this discussion that a contribution with AT = 2 in Eq. ilSSl) can only come from the alternative 2 in 
Eq. (lAllll . Therefore, setting L = L' = 1 and 1 = 1' = 0 in Eq. ifTmli . and making use of Eq. (6.3.2) of Ref. we 
get 


XQp = Ip,^; L = L' = 1) 
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(_y+z' 

—— 6vii 5u.,i' 

ZW 




A,7 A'.7' 


K 1/2 1/2 1 f K J’ J 
jp J J' J 15 A A' 


I V k ' 

1 1 K > 
X X K 

/ 


^ 1/2 1/2 1 f A ; 1 1 

> 1 J J J 5 j 

S' 1/2 1/2 1 f A' I' 1 1 

Jp 1 J' J \ J' 5 j ■ 


(A12) 


To proceed, it will be unavoidable to perform some manipulations with 12j symbols, and the needed identities are 
collected in Appendix ini for convenience. With the help of well known symmetry properties of 6j and 9j symbols 
|25|. one can make use of Eqs. inn), ins and ins), in succession, to perform, first the summation over A, and then 
that over X', in Eq. ISIS, getting 


XQp = Ip,^; L = L' = l) 


-- Ol^P Ol, ,1! 

311' 


k S S' \ 
S I' I j 


X 


,7 J' 


k S S'' 
K J J' > 
K 1 1 

/ 


K 1/2 1/2 

Jp J J' 
S' 1/2 1/2 
Jp 1 J' 


S 1/2 1/2 1 
Jp 1 J ] 


(A13) 


Repeating the same procedure, we can now perform, first the summation over J, and then that over J', to get 


XQp = Ip,^; L = L' = l) 


(_)^+fc+i'+S' 

Zip 


5uii Supi’ 


fc ^ 5' 1 

S I' I j 


X 


{ K 1/2 1/2 ' 

K 1/2 1/2 > 

k S S' ^ 


K 1/2 1/2 1 
Jp 1 1 J 


(A14) 


The 9j, as well as the last 6/, in Eq. IIA14II restricts AT to 0 and 1, and we conclude that the contribution with K = 2 
in Eq. (PI vanishes identically, both for k = 0, and for k = 1. Notice that this result holds, not only for jp = 1^3/2; 
which is of direct interest for but also for Jp = lpi/ 2 , which may be relevant for other p-shell hypernuclei. 


APPENDIX B: A = 1 CONTRIBUTION TO ljo 

Recalling Eq. (P| . it is clear that the A = 1 contribution to loq in Eq. 13HI, from the single-proton state Jp, has 
the form 


t^oGp, A = 1) 


= d cos 0. 


Pi 



pj dpi S 


P2_ , |Pl +P2p 
V2M 2M 2Mf 



X f{p,P)[Yi{9p,n)<»Yi{9p,0)]oo, 


(Bl) 
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where f{p,P) represents the rest of the integrand in Eq. IdSII . the important point being that it depends on the 
momenta only through p and P. 

From the explicit expressions of the spherical harmonics, we find 

/o 

[Fi( 6 >p, 7 r) 0 yL(6>p,0)]oo = - — cos( 6 »p + 9p), (B2) 

47r 

and, making use of the last two equations in ca, this becomes 

[n(0p,7r)®ri(0p,O)]oo = -^^|^. (B3) 

Introducing this result in Eq. m, we are left with 


‘^oOp, K =1) 


47r 


d cos 0. 


Pi 


X 


f{p,p) 


2 2 
P2-Pi 

2pP 



/ 


p\ dpi 6 



pI , \Pl+P2? _ ^ 
2M 2Mp ^ 


(B4) 


Let us now perform the interchange of dummy variables pi <-> p 2 in Eq. (IB4II . keeping cos Op-^ fixed. Then, noticing 
that, according to the first two equations in dia, p and P are invariant under this transformation, we arrive at the 
result 


‘^oOp, K =1) = - wo(jp, K = 1), 


(B5) 


from which it follows that the contribution with if = 1 in Eq. ilSSl) vanishes for k = 0. 


APPENDIX C: if = 1 CONTRIBUTION TO wi 


TABLE III: Results for the if = 1 contribution to oJi in the nonmesonic decay of ^aC. See text for detailed explanation. 


Model / Calculations 

t^l(lP3/2, K =1) 

AflA/ciA (%) 

TT 

Strictly local 

0.0007 

0.18 

Plus corrections 

0.0010 

0.20 

{■k,p,K) 



Strictly local 

0.0005 

0.15 

Plus corrections 

0.0007 

0.16 

TT + p 



Strictly local 

0.0007 

0.34 

Plus corrections 

0.0008 

0.26 

{n,ri,K) + {p,u},K*) 



Strictly local 

0.0003 

0.06 

Plus corrections 

0.0006 

0.12 


We have not been able to find an analytical proof that the if = 1 contribution in Eq. vanishes also for At = 1 . 
However, for the cases we are dealing with, this contribution can only arise from the jp = IP 3/2 proton state in 
and we have numerically computed its values in the several OME models we are considering. They are given, in 
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units of in Table HTTI The nonzero values may be due to truncation and roundoff errors. For instance, we have 
computed, with the same routine, the analogous contribution with k = 0 in the case of the complete OME model 
plus kinematical and nonlocality terms. Even though it has been proved in Appendix ^ that this is exactly zero, the 
numerical result came out as 0.0004, which is comparable to the value obtained for tt’i(lp 3 / 2 , AT = 1) in the same 
model. 

Comparison of Tables QH and cni immediately shows that the K = 1 contribution to oji is, in any case, very small. 
Furthermore, in the last column of Table ITTTI we give the relative effect that its inclusion would have on oa, and it 
always stays below 0.4%. We conclude that, even if this contribution is not exactly zero, it can be safely neglected. 


APPENDIX D: SOME PROPERTIES OF 12j SYMBOLS 


The 12j symbols arise in the recoupling of five angular momenta |27l l2Sj . They are not unique, but here we shall 
need only those of the first kind. 


Jl J2 J3 J4 

h h h Ia 
ki k 2 fca ki 


E(-) 


R4-X ^2 


jl ki X 

k 2 J2 h 


j2 k2 X 

ks J3 h 


j3 k3 X 

ki ji 13 


ji ki X 
jl ki I 4 


(Dl) 


as defined in Eq. (19.1) of Ref. whose notation for the 12j symbols we follow. In Eq. 1D1 |l . R 4 stands for the 
sum of all the angular momentum labels in the 12j symbol. 

These symbols obey the recursion relation Eq.(A.6.13)] 


Jl J2 J3 Ji 

h I 2 I 3 h 

ki k2 ks ki 


, _ ^_y2+k2+j4 + k4 



k2 ji X 
li li ki 


ki J2 X 
h li jl 


and have several symmetry properties, among which |28L Eq.(17.4)] 


f Jl 

J2 

j3 

ji ] 


\ 

js 

ji 

ki 1 

) 

h 

I 2 I 3 

h 

= 


h 

h 

h h 

[fcl 

k2 

k3 

ki J 

1 1 

[k2 

k3 

ki 

n J 


There are also reduction formulas, such as [2^ Eq.(A.6.39)] 


(D2) 


(D3) 


E 


= i-y 


Jl J2 J3 

Ji 


li h 

I 3 

li 

X k2 fca 

ki 


i+h+Ja+ks-ki ) 

J2 ki 

js I 3 

ki 

ji 


h ks k 2 


li k2 X 
ji li ki 

J2 ki ki 
li h jl 


(D4) 
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